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Abstract 

Parametric scaling representations are obtained and studied for the asymptotic behavior of inter- 
facial tensions in the full neighborhood of a fluid (or Ising-type) critical endpoint, i.e., as a function 
both of temperature and of density/order parameter or chemical potential/ordering field. Accurate 
nonclassical critical exponents and reliable estimates for the universal amplitude ratios are included 
naturally on the basis of the "extended de Gennes-Fisher" local-functional theory. Serious defects 
in previous scaling treatments are rectified and complete wetting behavior is represented; however, 
quantitatively small, but unphysical residual nonanalyticities on the wetting side of the critical 
isotherm are smoothed out "manually." Comparisons with the limited available observations are 
presented elsewhere but the theory invites new, searching experiments and simulations, e.g., for 
the vapor-liquid interfacial tension on the two sides of the critical endpoint isotherm for which an 
amplitude ratio —3.25 ± 0.05 is predicted. 

PACS numbers: 68.05.-n, 64.60.Fr, 64.70.Fx, 05.70.Jk 
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I. INTRODUCTION AND SCALING THEORY 



Consider, for concreteness, a binary liquid mixture consisting of two species, A and B. 
For a full thermodynamic description, one needs three field variables, say (T, h, g) ^,0]. The 
ordering field h is conjugate to the order parameter, M. For fluids, the order parameter may 
(in leading order) be taken as the number density Ap = p — pa(T, g) measured relative to a 
coexistence value po(T,g). Alternatively, M could be a composition variable such as mole 
fraction difference, a volume fraction difference, and so forth. For the nonordering field g, 
one may take the pressure, or the chemical potential of one species, either A or B, etc. 

Fig. ^a) illustrates a typical phase diagram in the three-dimensional field space Q. At 
any point on the surface labeled h = 0, the system exhibits phase separation into two 
coexisting phases, (3 and 7, rich in A and B, respectively. We will adopt the convention 
that the 7 phase has the higher (mass) density and hence sits at the bottom of a container 
when a gravitational field is present: see the inset in Fig. Q^a). By increasing temperature 
while keeping h = h — ho(T, g) = 0, the state point will reach the line A which is a locus of 
critical points, T c (g). Further temperature increase results in mixing of the (3 and 7 phases 
into a single phase, say /3j. On the other hand, decreasing g at fixed T < T c (g) on the h = 
surface leads to a triple-point line, r, at which appears a new, noncritical or 'spectator' 
phase a which represents the common vapor of the liquid phases (3, 7, and ^7: see the 
inset in Fig. ^a). A first-order transition, between the vapor and the liquid phases, occurs 
across the vapor-pressure surface labeled a which meets the h = surface at the triple line. 
The critical line, A, and the triple point line, r, terminate at a point (T e ,0,g e ): that is the 
"critical endpoint." 

Recent field-theoretic renormalization group theory has confirmed explicitly that the 
critical behavior at a critical endpoint is the same as on the critical locus 0,0]. Nevertheless, 
further, new bulk thermodynamic singularities do appear at a critical endpoint j^l 0, Q] . 

Beyond the bulk, however, there are singularities in interfacial or surface tensions when, 
in the presence of the vapor a, the two phases (3 and 7 merge into the homogeneous phase /3j, 
or vice versa 0, 0| . In Fig. ^a) one may follow the triple line and its smooth extension on 
the a surface, beyond the critical endpoint, or, in Fig. Hfb), simply trace the line g a (T). The 
"critical surface tension" between the coexisting phases vanishes below the critical endpoint 
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FIG. 1: (a) Phase diagram of a binary liquid mixture in the three-dimensional field space (T, h, g): 
see text for details, (b) Section of the phase diagram (a) containing the plane h = 0. The critical 
endpoint, (T e ,g e ), is where the critical line A terminates at the first order transition line g a (T). 

temperature T c = T e , as 

E 0h (T) « K\t\», t = (T-T e )/T e ^0-, (h = 0), (1.1) 

where, via standard scaling relations S, 0, 0, 10 1, the critical exponent is given by fi = 
2 — a — v so that /i ~ 1.26 for typical, three-dimensional fluids. The amplitude K has 
dimensions of energy per unit area where, here and below, we adhere to the notation set 
out in the Appendix of Ref. 8. The "noncritical tensions," £ a |/3 7 and Ti a \p, should behave 
, after subtraction of a suitable, nonsingular common background, E (T), as 



A£ a |/3 7 



K~\t\ 



0- 



Ch = 0-), 



[1.2) 

;i.3) 



which relations serve to define the amplitudes K + and K~ . Beneath T c the 7 phase coexists 
with the a and (3 phases and, in a sealed container, it sits below the f3 phase owing to its 
presumed heavier density. Now consider, hypothetically, bringing into contact the two phases 
a and 7, near (3^ criticality; this produces a new interface. The corresponding noncritical 
surface tension £ a | 7 can be obtained from Antonow's rule |l| which states 
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This relation can be derived by supposing that the three phases a, /3, and 7 can coexist 
and meet with nonzero contact angles, and then by letting the contact angle between the 
interfaces a\f3 and /3|7 go to zero (It should be recalled, however, that Antonow's rule 
typically fails at lower temperatures, specifically below a wetting temperature T w .) 



Using the classical van der Waals or Cahn-Hilliard theory and a model free energy 
of the Landau-expansion type, Widom [3| has studied various properties of the noncritical 
interfaces, such as that between a and (3, near the critical endpoint. Later, nonclassical 
critical exponents were embodied into the local free energy expression via postulated scal- 
ing forms Q]. However, the original theory of Widom and Ramos-Gomez led to an 
unexpected type of correction in the surface tension, namely a |t| 7 term: this is more singu- 
lar than the nonanalytic leading term whenever the spatial dimensionality, d, exceeds 
3 — fj P, |^ which cannot be considered acceptable. 

Fisher and Upton [5( pointed out that, near the critical endpoint, the amplitude ratios 

P=(K + + K-)/K, Q = K + /K~, (1.5) 

should be universal. They reported mean-field calculations 0, 0, [h]] yielding 

p = -|(V2- 1) = -0.20710... , Q = -v / 2, (1.6) 



which should be valid for d > 4. However, to obtain more realistic values for d = 3, Fisher 
and Up t0n g B panted p r e. lml „ ary cations using a, extend de Gennes-Fisher 
(EdGF) local functional theory |6j for fluid interfaces combined with a simple "interpolated 
linear model" for the equation of state (as described in Ref. 9). This approach provided the 
significantly different estimates 

P~0.1 2 , Q--0.83. (1.7) 

More recently, the EdGF theory has also been applied to critical adsorption problems [lsI ]. 

Our aim here, apart from estimating these universal ratios more precisely, is to calculate 
the noncritical surface tension in nonzero ordering field, i.e., AE a |^ 7 , AS a |^, and AS a | 7 as 
a function of t and h (or the order parameter, M), on the a surface in the full vicinity 
of a critical endpoint. This is a basic problem of interfacial thermodynamics first broached 



experimentally in pioneering work by Nagarajan, Webb and Widom 



14j . In fact, there is just 



a single function of two variables, AS(t, h), that is to be sought once a suitable background 



Eq(£, h) is subtracted from the total surface tension, say £ a u(t, h). Furthermore, in the 
first instance our main concern must be with the singular, critical behavior which we may 
confidently expect to be described in scaling form so that AS/|t| M is related universally 
(for d < 4) to the scaled combination M/\t\@ or, equivalently, to h/\t\ A , where, in standard 
notation, A = (3 + 7 = [35. 

Three decades after the development of renormalization group theory one might expect 
this problem to be susceptible to such an approach. Unfortunately, however, the still remain- 
ing difficulties, both technical and conceptual, are profound despite the progress reported, 
for example, in the review articles by D. B. Abraham, H. W. Diehl, and D. Jasnow in 
Vol. 10 of the series Phase Transitions and Critical Phenomena edited by C. Domb and 
J. L. Lebowitz (Academic, 1986) and subsequent developments (some of which are refer- 
enced in further detail below). Accordingly we report here on calculations based on local 
functional theories going back historically to van der Waals analysis of the critical surface 
tension E/3| 7 (T). Specifically, we pick up and develop the proposals of Fisher and Upton [|| 
who advanced, in particular, the EdGF theory which can consistently embody the correct 
nonclassical critical point exponents, especially rj > 0. 

Such theories rely on the availability of an accurate description of the bulk thermody- 
namic properties. To that end we will make heavy use of the parametric formulation of 
scaling theory in the neighborhood of a bulk critical point as extended to represent the true 



correlation length, £,oc(TJi] 
local functional theories 



and so provide a basis for calculating interfacial tensions via 



P, slO]. For completeness and ease of reference we recall the basic 
parametric expressions here. As standard, one first has 

t = rk(6), h = r A l(6), M = r f3 m(6), (1.8) 

where k(8) is an even function of the "angular" variable 9 with ko = k(0) = 1 and k(±6 c ) = 
so that 9 = ±6 C corresponds to the critical isotherm, T = T c , while 1(6) and m(9) are odd, 
with 1(0) = m(0) = and l(±0 x ) = with m(6A > so that 9 = ±9 1 describes the 
coexistence surface h = beneath T c : see Fig. 2 of [9J. 

For general thermodynamic purposes, however, it proves more effective to avoid integrat- 
ing the equation of state to obtain the free energy; accordingly [9j, we opt to treat h [and 
1(6)] as derived from the singular part of the reduced Helmholtz free energy which may be 
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written in scaling form as 

A s (t, M) = r 2 ~ a n(6), (1.9) 

where, following [9j, we can write 

n(6) = r- 2+a sing S^j h(M'; T)rfM'| , (1.10) 

in which the operation sing{«} extracts only the leading singular part while Mr < is a 
fixed reference value: see also ()2.1|) below. One then finds that 1(9) is readily expressed in 
terms of k(9) and n{6): see Eq. (4.4) of Q- 

The generalized local functional theories proposed in [f| also require the true correlation 
length, £,oo{T, h), which specifies the exponential decay of correlations (in the presumed 
absence of long-range power- law or van-der-Waals-type interactions). The corresponding 
scaling form can be written [l| 

ej2x = r-^ aoo (6), (1.11) 

where \ = (dM/dh)r is the reduced compressibility (or susceptibility) while r\ and v are 
the standard correlation critical exponents. It is worth stressing that an essential feature of 
the generalized local functional theories f| is to provide in a consistent way for r\ > since 
this is vital for the accurate description of real and realistic model systems when d < 4. 

The local functional theories, if they are to yield computations for the interfacial tension, 
also require that A S (M, T) and ^/2x are well defined through the two-phase region below 
T c where \M\ is less than M (T) « -B|t| , the order parameter at coexistence (or, in magnetic 
terms, the spontaneous magnetization). While this is most certainly questionable from a 
rigorous viewpoint, one may in practice construct trigonometric forms for k(9), m{9), n{9), 
etc., which extrapolate smoothly (and, indeed, analytically) to \9\ > Q\ and so through the 
two-phase region: see 0,0,0]. In such cases we take k(9), m(9), etc., as smooth periodic 
functions, of appropriate parity, in the interval —9q < 9 < 9q where #o then corresponds to 
h = M = for T < T c . On this framework an 'extended sine model' has been built and fitted 
to reliable estimates of critical exponents and amplitude ratios for the (d = 3)-dimensional 
Ising model 0\. The resulting scaling functions will be used here to study the interfacial 
tensions near a critical endpoint. 

In the scaling region the singular part of the full interfacial tension can consequently be 
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written parametrically as 

AE(t,/i) =7*3(6), (1.12) 

and our basic task is to calculate the angular surface tension function s(6). Note that AS 
represents (i) AS a m 7 when \9\ < 9 C , (ii) AS^ when 9 C < 9 < 9i, and (hi) AS Q | 7 when 
— 9\ < 9 < —9 C , in accordance with Fig. ^ and the notation explained above. Once s(9) is 
determined, the surface tension can also be written in the standard scaling form 

AS ps K\t\ ti S±(Dh/\t\ A ), (1.13) 

where the universal scaling function S±(x) can be readily calculated. (Note that, as cus- 
tomary, the subscripts + and — stand for t > and t < 0, respectively.) The amplitudes K 
and D here are the nonuniversal metric factors needed for normalization and to make S±(x) 
and the argument x, dimensionless. In the notation of [s| we take D = C + / B and x = h, 
which provided a convenient normalized variable in the analysis of Ref. 9. 



Experimentally, as mentioned, Nagarajan, Webb, and Widom (NWW) jlj] were the first 
to test theoretical predictions for the universal surface-tension scaling functions off the h = 
axis in their studies of mixtures of isobutyric acid and water. For the same mixture, Howland, 
Wang and Knobler [3| measured the critical surface tension S^ and Greer jl(| measured 
densities on the coexistence curve. The two latter experiments can be used to provide a 
consistency check and calibration of the NWW data |l7(. Other mixtures have also been 
examined. Quasi-binary mixtures of n-octadecane and n-nonadecane in ethane have been 
studied by Pegg et al. jl8( to measure the surface tensions through and near both the upper 
and the lower critical endpoints. The surface tension of the water and 2,5-lutidine system 
at and off the critical composition has been measured by Amara et al. |19J. For a similar 
mixture of water and 2,6-lutidine, Mainzer-Althof and Woermann determined the values of 
P and Q expe—* Q. Interfaeial tensions of the epical nnxtnre of 2-bnfoxyetnano. 
and water have been measured by Ataiyan and Woermann j^]. The applications of the 
present theory to these various data will be presented elsewhere j^. 

The rest of this article proceeds as follows. In Sec. II, we review briefly the classical 
theory of interfaces. The construction of more general local free energy functionals is taken 
up in Sec. III. Following Fisher and Upton we work out the details of the extended 

de Gennes-Fisher (EdGF) ansatz and obtain formulae for the equilibrium order parameter 



profile and the surface tension. The hypothesis that the noncritical vapor phase a can be 



replaced by a wall with a surface field hi p, |6| is introduced in Sec. IV. In the scaling limit 
hi/\t\ Al — ► — oo (where Ai is the appropriate surface critical exponent), there appear terms 
in the total wall tension that diverge although remaining analytic in t. After subtracting 
these divergent terms, we can express the finite singular part of the surface tension near 
a critical endpoint explicitly in parametric scaling form. In Sec. V these expressions are 
evaluated numerically for d = 3. However, an unphysical although quite small cusp is 
uncovered in the basic scaling function s(8) in (jl.l2j) . Its origin is discussed and found to 
reside in a fairly subtle deficiency of the EdGF scheme. By using an interpolation scheme, 
the cusp can be smoothed out leading to acceptable approximations for the universal scaling 
functions S±(x) in (|1.13|) . On this basis various concrete numerical results are presented in 
Sec. VI for the surface tensions in the vicinity of a critical endpoint. Sec. VI contains some 
brief concluding remarks. 



II. LOCAL FUNCTIONAL THEORIES FOR FLUID INTERFACES 
A. The Auxiliary Free Energy Function 

Let us consider various free energies that will be needed in discussing the local functional 
theory of fluid interfaces in a general way. Let A(M,T) be the true equilibrium Helmholtz 
free energy density that preserves the appropriate convexity properties in M and T j^. The 
free energy density A(M,T) = A tota i(M,T)/V can be obtained by integrating the equation 
of state, 

pM 

A(M,T)= h(M')dM' + A{M R ,T), (2.1) 

J M R 

using a fixed reference value Mr ^ 0: such a choice of reference value guarantees no sin- 
gularity across t = in A(Mr,T). In order to keep track of dimensions, we take M as the 
density difference (p — p c ) from now on. 

Now let A\M,T) be the Helmholtz free energy density in the one-phase region and its 



presumed analytic continuation into the multiphase region (f| ^, 0]. The typical van der 
Waals loop, which does not respect the convexity, should appear in the multiphase region of 
A'(M,T). The Maxwell construction applied to A'(M, T) repairs the convexity (although 
this is, of course, ad hoc). Evidently, A(M,T) = A\M,T) outside the multiphase region 
while A\M,T) > A(M, T) inside. The excess free energy A^(M,T) - A(M, T) which is 



S 



then always nonnegative, serves in the local functional theories, to determine the structures 



24] 



of the interfaces between the coexisting phases 

The conjugate free energy density F(hoo,T) = F tota i(hoo,T)/V, where the subscript oo 
denotes a bulk equilibrium quantity, may be obtained from A(M, T) via the Legendre trans- 
form, 

Fih^; T) = mm[A(M, T) - h^M], (2.2) 
= .1U/ X .V) -// x .\/ x . (2.3) 

where = M(T, h^) is the bulk equilibrium value of the order parameter. 

The interfacial tension that we aim to calculate is the excess free energy of a system in 
equilibrium created by one or more interfaces. All local functional theories for interfaces 
introduce an auxiliary free energy W[M(z);T,h oc ] that resembles the excess free energy: 
W(M) is always nonnegative and vanishes only when the profile, M(z), takes an equilibrium 
value of the order parameter corresponding to one of the coexisting phases (M Q for the a 
phase, Mp for j3, etc.). Since the equilibrium values of the order parameter M a , Mp, etc., 
are specified both by T and hoo, the dependence of the auxiliary free energy on T and hoo 
must not be overlooked. 

Thus the auxiliary free energy function W(M; T, hoo) needed in a local functional theory 



can be defined by 



15, 61 



W(M; T, h^) = A\M, T) - h^M - F{h^ T), (2.4) 
= A\M } T) - A\M X , T) - /ioo(M - MJ, (2.5) 

where we have used A* (Moo) — A(M 00 ). Since A'(M) represents the analytic continuation 
of A(M), we may also use ({23$ to write 

/•M 

W(M;T,h 00 )= / h(M',T)dM'-hoo(M-Moo), (2.6) 

J Moo 

because the equation of state must then equally have an analytic continuation: notice that 
A(Mr,T) cancels out. From this one immediately finds that W(M OQ \T ) h QO ) = and 
(dW/dM) M =M oa = 0. Starting from A^{M,T) > A(M, T) and using (Q and Q one 
obtains 

W(M; T, Zioo) > A(M, T) - hooM - min[A(M, T) - h^M] > 0. (2.7) 
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Thus, W(M; T, hoc) vanishes at the equilibrium values of M for all T and and its leading 
term in the expansion about is quadratic in (M— Moo). Otherwise, W takes only positive 
values. 



B. van der Waals Theory 

It is helpful to review briefly the van der Waals theory of interfaces which is a Landau-type 



classical theory 
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Assuming the existence of the local free energy [M(z); T, h^] 
that can be expanded in powers of M and T — T c , the van der Waals theory takes the local 
excess free-energy density functional as a sum of two terms [lj], namely, 

Af[M(z)} = W[M(z)] + | Jo (^fj 2 , (2.8) 

where Jo (= £ 2 /x) is a constant and z is the perpendicular distance from the interface pre- 
sumed to be flat. The square-gradient term, (dM/dz) 2 , accounts for spatial inhomogeniety 
in the simplest manner. The overall excess free energy is then given by a volume integral 
of Af(z). Translational invariance parallel to the interface enables one to factor out the 
area in the volume integral so that the free energy per unit area, or surface tension, can be 
written as jl| 

E[M(z)]= j ' dzAf[M(z)}. (2.9) 

Functional minimization of S[M] with respect to M(z) yields a differential equation for the 
equilibrium order parameter profile M(z). To supply boundary conditions let us consider, 
for example, a system containing two bulk equilibrium phases with M = Mp (or M_ OQ ) 
and M 7 (or M +00 ) located at z — — oo and +oo, respectively. For convenience, Mp and 
M 7 can be taken equal in magnitude but opposite in sign (M 7 > 0). Then, near criticality 
the equilibrium order parameter profile behaves like M 7 tanh(z/^) for the stated boundary 
conditions and the resulting surface tension is 

E= / dM ^2J W(M). (2.10) 
Jm s 

With a suitable representation for W[M], one can also study the critical and noncritical 



261 ] . However, this classical 



surface tensions Q,Q]> as well as the critical wetting transition 
square-gradient theory cannot satisfactorily embody all the correct critical exponents. This 
is because the square-gradient term in (|2.8j) . when the analysis is adapted to study the 

10 



decay of correlations (rather than just the overall interfacial free energy), always implies the 
exponent value rj — that is valid only for d > 4. Thus even upon using an equation of 
state that obeys scaling and embodies correct exponents (such as the extended sine model 
discussed in j2j), a square-gradient form for the local free-energy yields the classical value 
77 = for the correlation function decay 27 1. 

III. FISHER-UPTON THEORY FOR FLUID INTERFACES 

In order to generalize the classical square-gradient theory, we start with the total free 
energy of the binary mixture in the following form 

(T, hoo, g) = min F[M], (3.1) 
m (r) 

where T[M\ T, hoo] is a sought-for local free energy functional. Assuming translational in- 
variance along directions parallel to the interfaces, we may simply take M(r) = M(z). The 
solution M(z) that minimizes the functional T\M\ is the equilibrium order parameter pro- 
file for the specified T, hoo and boundary conditions. We can suppose that JF[M] has a 
homogeneous part describing the uniform bulk phases and an inhomogeneous part AjF, so 
that 

F[M{z)\ = F bulk (T, 9) + AJF, (3.2) 

where F bulk (T, hoo,g) is the bulk free energy (and g is the third field: see Fig. [TJ. If there 
were no interfaces, we would have AjF = 0, and Ftotai = ^buik- Following Jfj, we now consider 
a general local free energy functional of the form 

AF = J drA(M,M;T,hoo,g), (3.3) 

where M = dM/dz. Without great loss of generality, we may write 

A(M, M] T, hoo, g) = W(M)[1 + J{M)G{hM)\, (3.4) 

where J and A are to be functions of M and T. The function Q(x) should be even since 
the sign of M cannot matter in the free energy. Now A should vanish when M = so that 
AJ 7 = 0; this implies Q(0) = 0. To determine other properties of Q{x) we must proceed 
further. 
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For semi-infinite cases where an external wall is located at z = 0, the postulate ()3.3|) 
must be modified by adding a boundary term [28 \ 

AJF[M] 



A 



dzA(M,M) + f 1 (M 1 ;h 1 ), 



(3.5) 



where A is the area of the interface. The surface field h\ acts only at z = and is coupled 
to M\ = M(z = 0), the order parameter on the boundary. 
Now let us minimize 



&F[M{z)\ 



dz[A(M } M)+S(z)f 1 (M)) 



(3.6) 



A jo 

with respect to M(z) where the boundary term has been absorbed into the integral. Then 



in the usual way, integrating by parts, one obtain 



6AF[M(z)] 
A 



dz 



OA d { OA 



dM dz 



.dM, 



SM+5(z) 



djj_ 
dM 



OA \ 
dM J 



5M 



OA 

dM 



5M 



(3.7) 



and hence finds 



OA d I OA 



dM dz 



.dM, 



OA df x 



J)M dM / ,_ n 



0, and 



OA 



dM 



(3.8) 
(3.9) 



However, since A(M, M) has no explicit z-dependence, we can integrate ()3.8|) to get the 
first-order differential equation 



• dA 

A(M, M) - M — - = C, 
dM 



(3.10) 



in which C is a constant. 

Now consider a fully infinite situation so that the lower limit in ()3.5|) becomes — oo and the 
surface term drops out. Functional minimization again yields ()3.8)1 . and thence ()3. !())) ■ while 
the first (wall) condition in ()3.9|) becomes simply the bulk condition (8A/dM) z= _ 00 = 0. 

In the semi-infinite situation, the form ()3.4|) leads to 

dfi 



W{M l )J{M l )K{M 1 )g\KM) 



=o dM x 



(3.11) 
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which is an equation determining M\. Similarly, the two bulk conditions yield 
dA 



dM 



= W(M ±00 )J(M ±00 )A(M ±00 )g'(AM ±00 ) = 0; (3.12) 

±00 

but since W(M± ao ) = these conditions should hold automatically provided that J(M ±QO ), 
A(M±oo), and Q'(hM± oa ) do not diverge. We will see below that these functions have nice 
behavior so that we may forget the bulk conditions. 

Far away from the wall or the interfaces, in the infinite and semi-infinite cases, one should 
see only the bulk behavior of the system. Hence, it is natural to require 

M(z) -> M ioo , M(z) -> 0, as \z\ -> 00. (3.13) 

Because A(M± OD ,0) = the constant C in the first integral ()3.1()|) must thus vanish. Then 
using the postulate (|3.4|) one is led to 

which represents a differential equation for the equilibrium profile M cq (z). 

In order to devise ap pro priate expressions for A(M, M) which generalize the earlier 
de Gennes-Fisher theory |29|, applicable only at T = T c , Fisher and Upton 6] introduced 
the following desirable physical features or desiderata: 

A. The correct, non-classical critical exponents should be embodied, both in d = 2 
and in d = 3. Conversely, a reduction to van der Waals or Landau theory should be implied 
whenever the classical critical exponents are assumed. 

B. Near the critical point, all thermodynamic functions should satisfy the scaling 
and analyticity requirements: in particular, singularities (or nonanalyticities) should appear 
only at bulk criticality, i.e., when t = h = M ro = 0. 

C. For the semi-infinite system, the critical-point decay of the profile should behave 
in accordance with scaling as 12! 



m 

M c {z) ~ z~ p/v '. (3.15) 



D. When two plates separated by a finite distance L are immersed in a critical system, 
the order parameter decay ()3.15|) should have a correction factor so that 



MJz) ~ z-W" 



l+j 2 (z/Lf + ■■■ , L^oo, (3.16) 
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where d* = (2 — a)/u and j 2 is some coefficient. This prediction of the de Gennes-Fisher 



theory |2£ 



has been verified by several analyses including exact Isin g m odel calculations for 



d — 2 p, |6[ and field-theoretic calculations in e = 4 — d dimensions [31 



32 



33]. 



E. Away from criticality, the order parameter should always decay exponentially 

AM(z) = M(z) - Moo ~ e~ z/ ^ z -> 00, (3.17) 

where, clearly, £(T, h) denotes the true correlation length, £oo, but for brevity we will neglect 
the subscript 00. (Note again that we are assuming the absence of power-law forces or the 
likelihood that, if present, they do not enter explicitly into the asymptotic scaling functions.) 

F. For a finite critical slab, in the same situation as D, the critical profile will exhibit 
a minimum for similar boundary conditions satisfying M(z = 0), M(z = L) > 0, or a zero 
for opposing boundary conditions M(z — 0) > and M(z = L) < |34J]. It is expected that 
this profile, M c (z), behaves analytically near z = zq as 

M c (z) = h(z-z ) x '[l + k 2 {z-z ) X2 + •••], (3.18) 

where F(i) Ai = for the similar case, or Ai = 1 for the opposing case. F(ii) for the further 
exponents, one should have A 2 = 2, A 4 = 4, etc. 

G. Away from criticality a square-gradient expansion in the local free energy func- 
tional is expected to be correct and so should be reproduced by a satisfactory theory. 

H. To describe adsorption on a wall at z — 0, the theory should be consistent with 
the thermodynamic relation 

AM w & = -(§),' < 3 ' 19 ' 

where £ is the wall (or surface) free energy. 

I. The order parameter profile M(z;T, h, g; L) should be analytic in all noncritical 
regions. (Compare with B above). 

The starting point of the Fisher- Upton theory is the de Gennes-Fisher (dGF) ansatz 
for T = T c which in the general expression ()3.4j) is given by 0, 0], 

J = const., A(M) = £(M)/M, Q{x) = |x| 2 "^, 

with 77 = 2r}/{d* + r}). (3.20) 
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This form satisfies A — D and F(i). However, the dGF theory applies only at criticality. The 
extension proposed by Fisher and Upton (EdGF) uses x(M, T) = (dM/dh)T and postulates 



J=l, A(M;T, hoo,g) = \J £ 2 (M, T)/2x(M, T)W(M; T, ft^). (3.21) 
In order to satisfy E, one finds that the condition 

0(1) = 1 with 0(x) = xQ\x) - 0(z), (3.22) 

must be satisfied. Also, in order to satisfy F, G{x) should behave as 

G(x) =G + G^x^l + hx~ T + l 2 x- 2r + •••], 

with T = 2p/(P + v). (3.23) 

Finally, for small x — > the validity of a gradient expansion, G, can be seen to require 

Q[x) = x 2 + G 2 x A + G 4 x 6 + ■■■ . (3.24) 

Using the EdGF postulate (j3.21|) and the condition ()3.22|) . the solution of the general 
profile equation ()3.14j) reduces to x — ±1, or 

M = ±1/A(M), (3.25) 

where the signs ± must be chosen appropriately. Hereafter, we take the + sign, for an 
increasing profile as z — > +oo. The wall free energy then follows from ()3.5|) as 

poo 

Z = AF[M e(l {z)}/A= [1+Q{1)] / dzW(M) + f 1 (M 1 ), (3.26) 

Jo 

where fl3.3j) . ()3.21|) . (|3.25j) and Q{— x) = G(x) have been used. Using ()3.25j) once again, one 
can rewrite the semi-infinite integral as 

E= [1+0(1)]/ dM W(M)A(M) + /i(Mi). (3.27) 

Now the thermodynamic consistency condition, H, leads to |6j 

1 + ^(1) = 2 or 0(1) = 1. (3.28) 

It is remarkable that both the profile and the wall or interfacial free energy, S(T, h), do not 
depend on the details of G(x) 6j. However, one must note that the EdGF ansatz fails to 
satisfy I in certain situations |6|. 
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In order to repair this last problem, a generalized de Gennes-Fisher ansatz (GdGF) was 
devised jfj. It satisfies all of the desiderata A-I. However, in contrast to the EdGF theory, the 
profile equation, the surface tension formulae, etc., now depend on Q{x) explicitly. In order to 
obtain quantitative results, one must devise a representation for Q (x) that reproduces (|3.23j) 
and an analogue of (|3.24j) . One can indeed achieve this; but the resulting calculations 
become considerably more complicated than those for the EdGF theory. Accordingly we 
have explored numerically only the EdGF formulation. 

It must be recognized, however, that both EdGF and GdGF theories ignore capillary- 
wave fluctuations of a free interface, which are important for d < 3 Also, as indicated 
above, the analytic continuations of W(M), £ 2 (M)/2x(M), etc., into the multi-phase region 
have no known meaning in the strict sense of rigorous statistical mechanics. However, we 
expect that both theories will produce reasonably reliable results when fitted to exact d = 2 
and d = 4 results and good estimates for various d = 3 parameters, since they embody many 
correct physical features. 

Indeed, Upton's e = 4 — d expansion results for the universal amplitude ratio Q = 
K + 1 ' K~ |3Jj| demonstrate this point quite well. Using the field-theoretic approach to surface 
critical phenomena, he obtained the exact e-expansion Q = -y/2 + 1.521 257 e + C(e 2 ). 
Then, using the EdGF theory and the linear parametric model (which is known to be exact 
to order e 2 Q, El), he found Q = -V2+ 1.522 96 2 e + C(e 2 ). The coefficients of e differ 
by only 0.1%. 



IV. EDGF EXPRESSIONS FOR SURFACE TENSION NEAR A CRITICAL END- 
POINT 

We now discuss in more explicit detail the application of the EdGF theory sketched above 
to the vicinity of a critical endpoint. Allowing for a boundary term, which will be discussed 
further below, ()3.27|) can be written as 

E(T,M)=2/ dM ^/W(M)e(M)/2 X (M) + fi(Mi), (4.1) 

J Mi 

where the arguments and T are understood and (|3.21|) and (|3.28j) have been used. 
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A. Scaling Forms 



To embody the appropriate nonclassical critical exponents and satisfy the desiderata 

A and B we should, clearly, adopt scaling forms for W(M;T, h^) and the combination 

£ 2 (M,T)/x(M,T). As discussed, the required expressions must continue analytically (or, at 

least, sufficiently smoothly) into the two-phase region \M\ < Mq(T). We should also recall 

the necessity for including the further field g and the lambda line T = T c (g): see Fig. 
□ 

Following |9|, we thus introduce the dimensionless asymptotic scaling variables 

m = M/B\if, h=[h-h (T,g)]/(B/C + )\tf + \ (4.2) 

where the reduced temperature deviation is now 0, 0| 

t=[T-T c (g)]/T e , (4.3) 

while, for convenience, the tilde on h now denotes the fully scaled field and, for brevity, the 
phase boundary term, h (T,g), will usually be neglected below. As usual [q] B = B(g) and 
C + = C + (g) are, respectively, the critical amplitudes of the spontaneous order, M (T,g) « 
B\t\P, and the susceptibility, Xo(T,g) = (dM/dh)x>T c ~ C + /t 7 , on the h = surface above 

Then, following [5J, |6(, we can write the scaling form 

e(M,T,g)/2 X (M,T,g) « \M\-^Z ± (m), (4.4) 

which is crucial in going beyond van der Waals or Landau square-gradient theory because it 
introduces the small but positive (for d < 4) exponent rj. As above, the subscripts + and — 
will always denote t > or < 0, respectively. However, it must be realized here and below 
that the scaling functions Z + and Z_ and, likewise, others are, in fact, representations of a 
single, generally analytic scaling function continuing smoothly through t < 0. Thus, more 
explicitly, to ensure the analyticity of £ 2 /2x across the surface T = T c (g) or t = (recall B) 
the scaling functions Z± must have large fh expansions of the form 



Z ± (m) = Z™ 



i+Y / z-(±\m\- i/ r 



n=l 



(4.5) 



When T — > T c (g)± the terms in the sum clearly generate only the integral powers t n as 
required by analyticity when 
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The correct analyticity is most conveniently incorporated by using the parametric rep- 
resentations of the scaling functions recalled in Sec. I. The only new features required near 
a critical endpoint are the replacement of t by t, as defined in (J4.3)) and allowance for the 
(smooth) dependence of the nonuniversal factors = m(0) and Iq = 1(0) on the field g: see 
(II. 8|) and 38(. With this understanding the parametric forms (|1.8|) - (|1.1U|) will be adopted. 
Then one can express the coefficients entering ()4.5|) as 

^o°° = aoo(0c)[m(0 c )]^ = a ra (m c )^, (4.6) 
Z oo = [mjB] 1 ^ fa^ | rjum' c \ ^ ^ 




KAggf j 1 2(1 + w) ^ c < w^-fi-w) 

2 " 2[k' c f \ (3 aooc m c ^ 

I Q ooc _ a ooc fee + W ">J>^_K_ _ — c j 1 ^ 

where the prime denotes differentiation and for brevity we have used a' ooc = a^ )0 (^ c ), etc., 
and so on. 

The required scaling form for W(M; T, h^) must be somewhat more elaborate because 
of the additional dependence on h^: see the original definition ()2.5|) . (Indeed, the result 
presented in p is somewhat misleading since the dependence on was suppressed and the 
expressions given apply only for = 0.) Let us consider, first, the parametric represen- 
tation for W following from (J2.6)) with the aid of (jl.K)j) . Note that specification of T and 
hoc implies, via (jl.8j) . parametric coordinates and 8^ (for the corresponding bulk phase) 
while the variation of M at constant T and can be described by coordinates r and 6. 
Thus from (|2.6|) and ()1.8j) we obtain 

W{r, 9; roo , On) « r 2 ' a n(6) - r^n^oo) - r^/(^) [r^m(0) - r^m(^)]. (4.9) 

When t 7^ 0, the relation t = rk(9) = r OQ k(9 OQ ) leads to the desired form 

W(M- T, h^) « r 2 -<M0; 0^), (4.10) 

where, with /i^ = r^ ) /(6' 00 ) and f)1.8j) for M, we have 

-[n(0oo) - /(#ocMM] |*(0)/A;(0oo)| 2 ~ a , (4.11) 
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while the singular part of the Helmholtz free energy extended into the two-phase region, as 
discussed above, is given by (Hi) H- 

When T = T c we have 9 = 9^ = 9 C and w(9; 9^) simplifies to yield 



w c = w(9 c ;9 c 



(3mJ c /(2 - a), 
w' c = {dw/d9) e=eoa= e c = [(3m c l' c - (1 - (3) l c m' c }/{\ - a). 



(4.12) 
(4.13) 



Consequently w{9' ) 9 0O ) is continuous and smooth through T = T c (g) (or 9 = 9^ = 9 C ). 
However, the presence of the powers of A and (2 — a) in (j4.11|) shows that the curvature 
(d 2 w/d9 2 ) is not continuous through T = T c (g). Nevertheless, this lack of analyticity of 
w(6; 6*00) is not, of itself, expected to lead to corresponding nonanalytic behavior in surface 
free energies, etc., since the underlying bulk free energies and correlation lengths do vary 
analytically. 

Now we can express W in the alternative scaling form 



W(M; T, hoo) « \M\ s+1 Y±(m;h 



where, with A 00 ^) = -n(9 00 )/\k(9 00 )\ 2 ' a , we have 



(4.14) 



Y±(m; h c 



A^ih^lBml- 5 - 1 - h^M - Moo(T, h^M] 



-5-1 



+Y n c 



n=l 



(4.15) 



[This reduces to the Fisher-Upton expressions when h^ — > for t > and t < although 
Eq. (7) of 0, should be corrected by changing (±\y\)^ n ^ to read (±|y|~ 1 / /3 ) n .] On rear- 
ranging for \m\ — > 00 we have 



1 ± Yf°\fh\- 1/I3 - sgn(M)y^°|m|- A//3 + Y™ a \m\ 
+Y 2 QO \m\- 2/ P ± y 3 °°|m|~ 3//3 + • • • 



-(2-a)//3 



(4.16) 
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where the various coefficients are given explicitly by 



where m" 



Y%° 



1 2-c 



yoo 
J 



yoo 
1 1 



A^ihoz) + (5 2 /C + )/i 00 m 00 (/ ioo ) /B d+ X 



(5+1 v" 00 
> 



n c /(m c 
[m c /B] 



\5+l 



VP 



K 

[mjBfP 



— - (5 + 1) — 



2(1 - a) < < Q3- a)(l + 5) 







rir.m. 



n 



m' c k" 




m"(6 c ), etc. 



(4.17) 
(4.18) 
(4.19) 
(4.20) 



(4.21) 



B. Representation of the Noncritical Phase 

The scaling forms just discussed provide a satisfactory representation of the phases /3, 7 
and /?7 (see Fig. ^) in the vicinity of the critical line and critical endpoint, but they seem to 



give no account at al 



Fisher and Upton P, (| advanced a hypothesis, f2, which asserts that as regards the singular 
contributions near a critical endpoint, the noncritical phase a (typically a vapor when (3 and 
7 are liquids) can be replaced by a rigid, inert wall, say uj, characterized only by a nonzero 
surface field hi favoring, say, the bulk critical phase (3. 

If one accepts the hypothesis f2, one must consider a wall at, say, z = with a corre- 
sponding wall free energy as introduced in (|H.5jl which, following a Landau approach (2^, 
may be expanded as 

/i(Mi) = -h x Mi + c 2 Ml + • • ■ . (4.22) 

The boundary condition ()3.11|) determining M\ can be rewritten for the EdGF theory us- 
ing (02J), and (ET23J) as 

If the surface field near bulk criticality scales as hi ~ |t| Al when hi — > 0, this relation 
leads to Ai = i (2 — a — 77^) = ji — f3. However, it is known that the critical exponents for 
surface quantities such as hi and Mi are characterized by exponents A 1; f3i, etc., that cannot 
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of the noncritical or spectator phase, a. To overcome this draw-back, 



be derived simply from the bulk exponents 40|, |4jJ, |42| . Hence the exponent relation implied 
by (|4.22|) for small hi is not, in fact, valid. However, this failure of the local functional 
theory is of no concern in the present situation since we wish to keep the field hi fixed and 
nonzero. Hence the correct scaling combination hi/\t\ 1 diverges when t — > and hence 
hi does not appear explicitly in the asymptotic scaling forms. From our viewpoint, then, 
it is only necessary to investigate large \hi\ and, correspondingly, large |Mi|. In a more 
general formulation the a phase would be represented by an extra ( "third" ) minimum in the 
extended free-energy function A^(M,t,g) at some value M = M a (t,g) with, say, M a < 
corresponding to a vapor phase: a fuller investigation then shows that the hypothesis f2 
is justified within EdGF theory. Of course, the effective surface order Mi does not become 
indefinitely large in this formulation: however, for the singular behavior near the critical 
endpoint we may imagine taking Mi to — oo, the negative sign being chosen so that, as 
indicated, the density p a of the vapor phase lies below the critical density p c (corresponding 
to M = 0) so that the wall field hi should induce a negative value of M. This limit entails, 
as will now be discussed, the subtraction of appropriate leading terms that would otherwise 
appear as divergencies but which, in reality, contribute only to noncritical 'background' 
terms. 



C. Implications of an Unbounded Surface Field 



To use the integral expression (|4.1|) for the surface tension in the limit where hi and, 
hence, M\ remain finite so that the scaled variable m ~ M/\tf becomes large near the wall 
at z = when t — > 0, we need to examine the integrand. When \m\ — > oo we find 

e {My 1/2 



W{M) 



2 X (M) 



\M\ 



G*-/9)//9 



/yoo 70c 



.1/2 



1 



-sgn(M) 



h^/BSy™ 



+ 



(Y{° + Z™)t 
2\M/B\ l /P 

v^oo U|2— a 
I 2-a\ L \ 



2\M/B\ S 2\M/B\ S+1 
[4(K>°° + Z 2 °°) - (Yf° - Z™) 2 ]t 2 



%\M/B\ 2 IP 

S 1 ' 4\M/B\VP+ S 
(Z™ - Y£°)Y™ a t\t\ 2 - a 



4|M/S| 1 //s+«+i 



+ 0(\m\~ 3 ^) 



(4.24) 
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On integration the leading term \My JJ '-p)IP yields | ik^! | which diverges when Mi — > — oo 
since fi//3 > 0. Similarly, the second term yields a divergence when (// — l)//3 > 0, which 
applies when d > 2. The third term gives |M 1 |^ _A )/ /3 which vanishes in the limit Mi — > — oo 
since /i < A for d = 3. However, when d = 4 it diverges as In \Mi\. Higher order terms do 
not yield any divergencies when Mi — > — oo. Hence, to remove the divergence in the integral 
when Mi — > — oo for 2 < d < 4, the first three terms must be subtracted. 

When Mqo > the integral in (J4.1j) runs through M = which is problematical because 
the piece to be subtracted is singular at M = 0. To avoid this, we split the integral at an 
arbitrary point M* < and make the subtraction only in the integral from Mi to M*\ of 
course, the values of M* should not affect the final results. Finally, the subtracted form of 
the surface tension expression (j4.1j) can be written 



E(t, /ioo; Mi) = h + I 2 + I 3 (M*) - h{Mi) + fi{Mi] 



(4.25) 



where the various contributions follow from 

-Moo 



h 
h 



2 / dM [W(M)e(M)/2 X (M)] 

t-M* 

2 / dM 



1/2 



Mi 



{W(M)e(M)/2 X (M)} 1/2 - {Y™Z™) 1/2 \M\ 



,1/2 I M \p/fi-i 



1 (YT + Zf)t 



2\M/B\yP 2\M/B\ l 



h{Mi) = 2(F °% 



Mi 



dM \M\ 



n/p-i 



(4.26) 



(4.27) 



(4.28) 



2\M/B\VP 2\M/B\ S 

The third, indefinite integral, which together with fi(Mi) contains all divergences in the 
limit Mi — > — oo, can be performed analytically yielding 



h{M) 



-2 (Y~Z£ 



I*/ p 



2y °° 



20* - 1)//? 
j3(—M)^ l3 ~ 5 /(fi — A), rf<4 
ln[(-M)/fl7£], d = 4 



(4.29) 



where for d = A we have made the argument of the logarithm dimensionless and scale- free; 
this is harmless because the extra ln(Sr^) term amounts merely to an additive constant. 
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D. Singular Part of the Surface Tension 



Following the discussion, we now identify the finite, singular part of the surface tension 



as 



AE(t,h c 



h + h 



with 



h{M*) 



(4.30) 



where the limit M\ — ► oo in any remaining Mi-dependence is to be understood. Note, 
furthermore that when d < 4, the diverging parts that have been subtracted vary only 
linearly with t. As mentioned, they may thus be regarded as a part of the common analytic 
background, S (T), for the surface tension: see (jl.2j) and (|1.3|) . In identifying AE for d = 4 
there is some unavoidable arbitrariness associated with the introduction of the ln(Sr^) 
term: but this is of little significance. 

Finally, to employ the parametric representations, we change the integration variable 
from M to 9 in the integrals I\ and l<i- From (jl.8|) . we have 



dM 



00 



and on defining 9* < via 



m fc(fl)m'(fl)-/?fr(fl)m(fl) u|g 



M*/\tf = m(9*)/\k(9*)f 



we find 



h = 2sgn(t)|t|^ 
I 2 = 2sgn(t)|t|" 



)Q k(9)m'(9)-Pk>(9)m(9) 
d9 rrr^ ^a 00 (9)w{9), 



d9 



\k{9)\^ 

k{9)m'{9) - j3k'(9)m{9) 
\k(9)\»+ l 



(4.31) 



(4.32) 



(4.33) 



y/a 00 {9)w{9) 



2\m{9)/B\ 1 /' 3 2\m(9)/B\ 5 
where w(9) is given by (j4.11|) (for fixed #oo), while J 3 (M*) can be written as 



(4.34) 



-2 (Y^Z, 



ooU/2 
J 



*\-fi 



/i 2{n-l)/P 



+ 



VOO 



/%*) A /Gu-A), 



d < 4 
d = 4 



(4.35) 



where y* = \k(9*)\[— m(9*)/ B]^ 1 ^ 13 . It is to be understood that, below T c , the integral I\ 
must be interpreted as = J e f° + J^ 00 since, by our periodic parametric construction, 
9 = +9q and 9 = —9 are identified together with M — 0. 



23 



These expressions for I l7 I 2 , and J3 together with ()4.3U|) for AT,(t, hoo) constitute our 
explicit results for the singular part of the surface tension. Note that Mi, the boundary value 
of the profile on the wall that stands in for the a phase, does not enter these expressions. 
On the other hand the arbitrary value 9* does appear; however, the total Ii+I 2 + 1% should 
be independent of 9* and this may be checked in explicit calculations. 



V. PARAMETRIC SCALING FUNCTION FOR d = 3 

In this section we examine the numerical results for the parametric scaling function s(9) 
introduced in (jl.l2|) that follow from the EdGF surface tension expressions obtained in the 
previous section, specifically 1)4.30)) with ()4.33|) — f)4.35j) . 



A. Numerical Evaluation of the Scaled Tension 



Now given the parametric angular functions k(9), m(9), etc., the numerical integration 
of I\ in ()4.33)) can be performed readily. However, there is an endpoint singularity in 1% 
in ()4.34)) at 9 = —9 C where k{9) vanishes. On using the expansion ()4.24|) for [W^ 2 /2x\ 1 ^ 2 
one sees that the singularity has the form \9 — 9 C \^ with 0=1 — a — \x = v — 1<0 for 
d = 3. Although divergent, this singularity is integrable and can be handled by standard 



techniques 



43]. 



As a check on the numerical calculations, it is useful to recall the mean-field limit for 
which exact analytic results can be obtained. Indeed, from the asymptotic classical equation 
of state, namely, 

h(M,t) = DM(B 2 t + M 2 ) with D = 1/B 2 C + , (5.1) 

one can find the auxiliary free energy via ()2.6|) to obtain 

W(M) = \tB 2 D{M 2 - M 2 J + \D(M 4 - M£) - h^M - M^), (5.2) 

in which = M(T, hoo) is, of course, the appropriate bulk equilibrium value. Using the 
EdGF formulation and taking £ 2 /2x = |J [see ()2.8|) ]. the critical surface tension can be 
calculated analytically, yielding 

£/3| 7 = 2 / dM [W(M)e/2x\ f (5.3) 
J-M Q (T) 



|V2Jo/5 2 C+M Q d (T), 
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from which the amplitude K defined in (jl.lj) is 

K = \K Q with K = y/2J Q B*/C+. (5.4) 

Similarly, the amplitudes for the noncritical surface tension can be calculated using the 
subtraction scheme discussed above; they are found to be 

K+ = -IV2K , K~ = \K Q , (5.5) 

which values confirm the universal amplitude ratios P and Q stated in ()1.6|1 . 

For the numerical results presented below we have used the extended sine model ex- 
pounded in [jj, specifically in Eqns. (5.6) and (7.4) with (6.1) and (4.4). The numerical 
parameter values for d = 3 are given in Eqs. (6.2) and (7.5) of Ref. 9. In light of more 
recent estimates for the critical exponents and universal amplitude ratios of d = 3 Ising mod- 
els 44,B Q, 

the model parameters should, ideally, be updated. However, the resulting 
changes could induce only rather small effects. Thus the new estimates v ~ 0.6302 and 
7 ~ 1.2375 45] differ from the adopted values of j^J by only 0.16% and 0.28%, respectively. 
The universal amplitude ratio estimates ft/fx = 1.963 ± 8 and C + /C~ = 4.762 ± 8 |45j 
agree well with those used in |9j lying within the uncertainty ranges although their cen- 
tral values are shifted by 0.15% and 3.8%, respectively. Furthermore, it has been shown 
via parameter-sensitivity checks that the associated variations in universal ratios are less 
than 5% 0. Thus the accepted parameter values from [l| are perfectly reasonable for the 
numerical aspects of the calculations reported here. 

The computation may be set up to generate the mean-field results applicable when d > 4: 
see Appendix D of Ref. 17 for the appropriate parameter values. Our numerical values for 
the mean-field case agree up to eight digits and the lack of any dependence on 9* was 
verified. (For convenience the values 9* = — 3# c /4 for < 9^ < 9 C and 9* = —(9 C + 9 )/2 for 
9 C < 9oo < 9 were adopted jl7j|.) In d = 3 dimensions the value of 9* was varied from —0.05 
to —0.26 for a sample value of 9^ near 9 C : no change in s(9) was found within eight-digit 
precision. 
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FIG. 2: The angular scaling function s(9) for the surface tension in d = 4 dimensions as given by 
EdGF theory. Some significant numerical values are 6 C ~ 0.562 345 and 6\ ~ 0.667 708 correspond- 
ing to s c /ao = 3.120 23, s- c /a = -0.122 322, Sl /a = 0.976 225, and s_i/a = 0.325 408. Note 
also s(0)/a = -3.306 20 |lj. 

B. Calculated Scaling Functions 

Our calculations may be summarized by presenting the angular function s(9) for the 
surface tension. For the purpose of normalization, we employ 

e>o = (floo^mo) 1/2 , (5.6) 

which has the same units as surface tension. Fig. El presents a plot of s(9) for d = 4. One 
sees that, owing to the wall dependence (embodied in Mi — > — oo), s(9) is not symmetric 
with respect to 9 = 0. Notice also that s(9) varies perfectly smoothly across the critical 
isotherm values 9 = ±9 C as it should. 

The qualitative features of the scaling function in d = 3 dimensions are the same as for 
d = 4. However, as seen in Fig. El an unexpected, albeit small cusp appears at 9 — 9 C . The 
presence of the cusp (which is certainly absent when d — 4) implies that there is a line of 
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nonanalyticity along the critical isotherm t — when > 0. This is quite unphysical since 
singularities in the thermodynamic functions, including surface quantities, can occur only at 
the critical point (t, h^) = (0,0). Thus, from the analyticity away from the critical point, 
one expects 

AE(t, /ioo) = AE^/ioo) + AEfih^t + AZfih^t 2 + ■ ■ • , (5.7) 

when /loo 7^ and t — > 0±, respectively. However, analysis shows that the surface tension 
predicted by the EdGF theory behaves as 

AS(t,/ioo) = A£ c (/ioo) + A£f(/ioo)|tr 1+a/2 

2 

+A£±(/ ioo )|tr 1 -"/ 2 - A + • • • , /ioo > 0, (5.8) 

when t —>■ 0±: see Appendix A in Ref. 17. In d > 4 one has \x — 1 + = | so one should, in 
principle, see a square root cusp then; but such a cusp is absent in Fig. El This is because the 
amplitudes AEf (/ioo) vanish identically in the classical situation: see Eq. (A. 30) of Ref. 17. 

2 

It transpires, as we will now explain, that this erroneous singular behavior appears be- 
cause the EdGF theory does not satisfy the desideratum I of Sec. III. More specifically, a 
zero of the order parameter profile is not represented correctly when T ~ T c and > 0: 
recall that the order parameter in the EdGF surface tension integral (|4.1|) now runs from 
M = Mi < to M = Moo > and so passes through the value M = which is critical 
when T = T c or 9 = 6 C . Equivalently, the EdGF surface tension can be represented by a 
real-space integral involving the free energy W[M(2)], as in (J3.26j) . Now, let us focus on 
the term h 00 M{z) in W(M) [see f)2.5|l ]: i.e., consider the contribution of jdzh QO M{z) to the 
surface tension. In order to find the scaling behavior of the profile near M = 0, we may 
analyze (|3.25|) to find (see Appendix B of Ref. 17 for details) 

M(z) (x-t nu/2 {z- z ), (5.9) 

when z —> Zq. Evidently the z-dependence is linear as expected from the analyticity of M{z) 
when t ^ 0. However, by scaling z varies as £ ~ t~ v and hco as t A ; thus the scaling behavior 
of the profile M(z) near z = z can be estimated leading to 

f dz h^Miz) ~ felloe f „ £ (W2)-2,+A_ (5 1Q) 



27 



2 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 

—0i — 6 C 6 9 C 6\ 



FIG. 3: The calculated angular function for the surface tension in d = 3 dimensions based on 
EdGF toy a„ d U- eXtend e d sin e m o d e> Q. T h e ^ sh ows a CUSp * . - (i.e. a t 

T = T c for /ioo > 0) with a corresponding value s c /o"o = 3.91731. The values of s_ c and s±i are 
entered in Table I. 

Note that the second integral is dimensionless and scale-free and so has been regarded 
as a fixed number. Via the scaling relations, the exponent {j]v/2) — 2v + A reduces to 
\i — 1 + (a/2) which corresponds to the cusp appearing in (|5.8|) . 

Since the correct scaling function should not exhibit such nonanalyticity, we introduce 
an interpolation procedure that smooths out the calculated cusp in s(8). The 9 range, say 
[6 a ,0b], selected for the interpolation is arbitrary. However, the choice (0 a ,0b) = (0.20,0.35) 
encompasses the cusp and seems reasonable. We have used a polynomial P(9) of degree 
5 which is the minimal order to match the first and second derivatives at both ends of 
this interval. The resulting polynomial for the extended sine model parameters is given 
in Eq. (4.4.50) of ^|. It agrees closely with the calculated values of s(9) in the intervals 
[0.20,0.26] and [0.32,0.35]; Furthermore, the largest deviation of P{9) from s{9) occurs at 
the cusp and is only 3%. 

From now on, when the distinction matters, we will write s(9) for the ameliorated angular 
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TABLE I: Numerical values for the ameliorated angular surface tension scaling function s(9) in 
three dimensions. Note that 9 C = 0.269 293 and 6>i = 0.422 519 j^. 



e 







s(0)/<ro 


e 


a(6)/<r 


9 


s(9)/<t 


-9! 


1.72713 


-0.20 


-1.67743 


0.08 


-0.07619 


0.30 


4.458 74 


-0.40 


1.125 36 


-0.16 


-1.783 25 


0.12 


0.58732 


0.32 


4.792 82 


-0.36 


0.20779 


-0.12 


-1.76707 


0.16 


1.355 90 


0.34 


5.022 93 


-0.32 


-0.513 65 


-0.08 


-1.638 86 


0.20 


2.21093 


0.36 


5.12125 


-0.28 


-1.055 67 


-0.04 


-1.405 58 


0.24 


3.129 71 


0.38 


5.049 63 


-e c 


-1.173 07 





-1.06944 


+e c 


3.812 64 


0.40 


4.739 23 


-0.24 


-1.438 71 


0.04 


-0.62787 


0.28 


4.05087 


+01 


3.806 99 



function to distinguish it from s(9) that has the cusp at 9 = 9 C . Selected numerical values 
of s(9 ) are given in Table I; values spaced at intervals A9 = 0.01 are available in Table 4.1 
of 

In the previous study of the extended sine model various parameter sets near the 
preferred set in Eq. (6.2) of [9J were examined in order to check the sensitivity of the univer- 
sal amplitude ratios. Here, we also check the sensitivity of s(9) under the variations in the 

n n 

parameter sets considered in Table I of [9j. As explained in Sec. VII of |9j, the optimal pa- 
rameters aooo, aoo2, etc., for the true correlation length are determined separately for each set 
by fitting the universal ratios K{f~) 2 and aA + {f + ) d , f + /f~, and (C + /C c ) {f c / f + f~ v . 
Then s{9) is computed for each parameter set using EdGF theory. 

It must be noted that the value of 9\ depends somewhat on the parameter set. This implies 
that the values of s(9) calculated from different parameter sets are not strictly comparable. 
However, the change in 9\ occurs only in the third decimal place. Thus, ignoring the small 
changes in the 9 scale, we have examined As (9) = s(9) — s(6%, i.e., the deviations of s(9) 
calculated using the sets 1, 2, . . . , 6 in Table I of |9[ in place of the preferred set (6.2) of 
Ref. 9. 

The largest deviation occurs at 9 = 9\ ~ 0.42; normalized by the value s(9\) in Table I 
it is about 3%. In |9j the variations in the predicted universal ratios associated with the 
parameter changes ranged from 0.2% to 5%. Thus, the function s(9) shows only the same 
level of variation as might have been anticipated. 
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Finally, we have also checked the effects on s(9) of using "untuned" estimates for 000(6*). 
Indeed it was found in 9] that the [2/0] Pade approximant for a OQ (9) could not fit the uni- 
versal ratio K(f-f while also fitting /+//", aA+(f+) 3 , and (C+ /C c )(f c /c + ) 2 ~ ri . Likewise 
for the other approximants of the same order, namely, [1/1] and [0/2]. Hence, it was neces- 
sary to introduce the [3/0] approximant. By using the [2/0] approximant for a OQ (9), which 
is representative of the low-order approximants, we have calculated the angular function 
s(9; [2/0]) with the preferred parameter set (6.2) of 9]. The difference between s(9; [2/0]) 
and s(9; [3/0]) in the range —6\ < 9 < 9 C is very small; indeed, it is invisible on a plot. How- 
ever, a large difference of about 30% arises in the subcritical range 9 C < 9 < 9\. This occurs 
(i) because the approximants [2/0] and [3/0] behave quite differently in the two-phase region 
9\ < \9\ < 9q (see Fig. 8 of [9() and (ii) because the calculation of sij)^) for 9 C < 9^ < 9\ 
involves integration through the two-phase region whereas the range —Q\ <9 00 <9 C requires 
integration only through the one-phase region in which the approximants are very similar. 

VI. NUMERICAL RESULTS FOR THE SURFACE TENSIONS 

Based on the angular scaling function s(9) or, where appropriate, the ameliorated version 
s(9), calculated with the preferred parameter set (6.2) and (7.5) of Ref. 9, we now describe 
various theoretical predictions for the interfacial tension near a critical endpoint. 

A. Amplitude Ratios 

Recall first the surface tension amplitudes and K defined in (|l.ljl - p.3|l . On using 
the parametric form for AS in ()1.12j) the surface tension amplitude above T c is readily read 
off as 

K + = s{0)/[k{0)]^ = s(0). (6.1) 

Owing to the asymmetry of s(9) with respect to 9 > 0, we must choose the sign of ±9i 
appropriately for zero-field below T c . In accord with Fig. GJa), the a\(3 interface exists 
beneath T c when h < 0; thus, the amplitude K~ must be evaluated at 9 = —9\ yielding 

K~ = s(-9 1 )/\k(9 1 )r. (6.2) 
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Similarly, at 9 = 9\ one obtains the amplitude for AE a | 7 in the limit ft — > 0+, which via 
Antonow's rule is the sum (K + K~): thus we have 

K=[s(9 1 )-s(-O 1 )]/\k(0 1 )\». (6.3) 

It is also instructive to define corresponding critical surface tension amplitudes on the 
critical isotherm. Via scaling one can write 

AS « Kl\hf A « KllMl^, (6.4) 

where the sub- and superscripts > and < stand for h, M > and ft,, M < 0, respectively: 
owing to the asymmetry of the surface tension with respect to ft — > —ft, this distinction is 
essential. For reference we may note that ///A ~ 0.806 while fi//3 = 3.85. The parametric 
representations then yield 

^> = s(^ c )/|Z(^)|^ A , if< = S (-^ c )/|/(^)r /A . (6.5) 

Note the appearance of s(9) here so that the prediction for K> depends on the amelioration 
procedure. 

The specific reduced surface tension amplitudes predicted by our EdGF theory with the 
extended sine model are thence 

K + /a = -1.069 44, R-/a = 1.281 68, K/a = 1.543 44, 

K >/ cr ol >l/A = °- 293 755 > K </^o l o >l/A = -0.090 382 6, 

K^/aorriQ^ 13 = 923.242, K^/a^" 3 = -284.063, (6.6) 

where we used Table I for the s(9) values and the extended sine model values [^] 

k(9x) = -1.266 16, l{9 c )/l = -l(-9 c )/l = 0.529 162, 

m(9 c )/m = -m{-9 c )/m = 0.240 366, (6.7) 

while o"o is defined in (15.6)1 and represents the nonuniversal surface tension scale. 
From Eq. ()6.6)1 . one obtains the universal amplitude ratios 

P = 0.137 5 ±2, Q = -0.834 ±2, (6.8) 

where the uncertainties have been estimated by examining results from the other parameter 
sets in Table I of Ref. 9. In comparison with the preliminary calculations quoted in (jl.7|) . 
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one sees that our improved estimate for P is about 15% larger [and still positive in con- 
trast to (jl.6|) ] while the Q value displays only a 0.5% deviation. These estimates will be 



22j in relation to the experimental observations of Mainzer-Althof and 



discussed elsewhere 
Woermann 

For the surface tensions on the critical isotherm, our analysis generates the universal 
amplitude ratio predictions 





( — 


K 






(± 


K 





3.34 ±4, -^LB^ 13 = 5.26 ±7, 
K 



| =-1.02 9 ±l, -J*i^ //3 = -1.62 ±2, (6.9) 

K 

where we recall that B and C + are defined as in |8j. More directly we find 

K>/K< = K^/K< = -3.25 ± 5, (6.10) 

which may also be used in analyzing experimental data and might well be tested in simula- 
tions of Ising-type systems. 

We plan, as mentioned in the Introduction, to consider the applications of the present 



theory to experiments in the future |22j; but it should be noted here that the difference 
in sign and magnitude of the amplitudes and was explicitly remarked by NWW 
on the basis of the theory of Ramos-Gomez and Widom [l^], a square-gradient approach 
formulated to incorporate 5 = 5. [See also Rowlinson and Widom [l] (pp. 287-293) and our 
comments following p. 41) and (2.10) above.] Furthermore, their analysis leads to K^/K^ ~ 
-148/42 ~ -3.5 2 : see |lo| p. 614 and |l| Eq. (9.124). This value is only some 8% larger in 
magnitude than found here. 



B. Scaling Functions 

While parametric scaling forms are conceptually and computationally effective, direct 
scaling representations, as in ()1.12j) . are more useful for comparison with existing or proposed 
observations. Thus from the angular function s(9) we have computed the scaling functions 
S±(h): in Fig. IHa) these are plotted in terms of the field (or chemical potential) variable h. 

In experiments, for example on binary mixtures, the density deviation (oc M) is more 
readily accessible than the chemical potential deviation (oc h). Hence, in practice scaling 
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FIG. 4: The universal scaling functions predicted by the EdGF surface tension theory (a) vs the 
ordering field h and (b) vs the order parameter in. Note the lack of symmetry about h = in = 
and the jumps dictated by (|6.12|) . 



33 



plots in terms of m are more convenient: the corresponding scaling functions may be defined 
via 

AE»if|f4(m) with m = M/B\t\ p , (6.11) 

while their behavior is shown in Fig. HJb). The normalizations adopted in ()1.12|) and here 
lead, with Antonow's rule (jl.4j) to the "jump conditions" 

S_(0+) " S-(O-) = S M (1) - S M (-1) = 1, (6.12) 

that enter below T c . 

Away from the critical point, all bulk and surface quantities must be analytic unless some 
phase boundary intervenes. This is reflected in the smooth behavior of S + (h) for all h. In 
contrast, there is a jump in S^(h) at h = in accordance with (|6.12j) : for h ^ 0, however, 
S-(h) is analytic for all h. Similar considerations apply to the break in the S^(rh) plot seen 
in Fig.0fb). On the critical isotherm T = T c the two branches of the scaling function, S + (h) 
and S-(h), must join smoothly when h — > ±oo; this again is a consequence of the overall 
requirement of analyticity away from the phase boundary. Hence, the two branches of S±(h) 
and S^{rh) must both asymptotically approach one another when h — > ±oo and m — > ±oo: 
this can be seen easily in Fig. HJb). Notice that the dotted line plots labelled T = T c in 
Fig. |3] depict only the asymptotic power-law behavior embodied in (|6.4j) that corresponds 
in these scaled plots to h — > ±oo. (It may also be remarked that the unphysical cusp in the 
EdGF surface tension on the critical isotherm for h > is located at h — > +oo; thus, with 
or without the cusp, the scaling plots in Fig. HJlook similar.) 

Finally, it should be clear that the scaling function plots in Fig. 0] may also be read as 
describing the variation of AS(T, h) oc S±(h) = S^irh) with h and M at fixed values of 
T < T c . In particular, one may then notice that the isotherms of AS vs M for T above T c 
will cross the critical isotherm when M is positive, and, by continuity, hence also cross one 
another This crossing of the surface tension isotherms above T c has been anticipated 
theoretically by Ramos-Gomez and Widom [10: see text after their Eq. (3.15) and Table I] 
who also point to some experimental evidence of crossings. 

On the other hand the isotherms of AS vs [h — h (T,g)] (which, recalling (|4.2j) et seq., 
should replace h in the figure) do not cross: see also Fig. 4.10 of [l2(. However, in addition 
to the displacement h (T, g), inclusion of the surface tension background S (t, h) will distort 
naive expectations based on Fig. 0] when real isotherms for the total interfacial tension are 
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examined vs density or chemical potential. See Fig. EJb) below for another aspect of this 
issue. 

It is appropriate to mention here that in their theoretical analysis of a binary fluid mix- 
ture, NWW [lj] tacitly assumed symmetry in the surface tension above and below T c by 
supposing S^(m) = S^(m): see their Eqs. (2.8) and (2.9) and Figs. 3 and 10. However, 
on physical grounds such a symmetry is quite implausible. Thus, below T c there are two 
distinct fluid phases (3 and 7 favored by h > and h < 0, respectively, and a vapor phase a 
(or wall) favoring the (3 phase; but, above T c , there is only one fluid phase (3^\ see Fig. [TJ 
Hence, one must allow for the T < T c symmetry breaking differences seen in Fig. 0] For 
comparison, one may recall that for the bulk equation of state when expressed in terms of m 
(or in terms of h) one also needs two scaling functions, say Q±(m) for t > and t < 0, as is 
well known; and, again, the two branches must be analytically related. However, in the bulk 
thermodynamics one may accept full asymptotic symmetry under h —h orm<^> — fh and 
then employ just a single function, say, in terms of the variable t/lml 1 ^ oc m~ 1///3 . But that 
symmetry is also not applicable to the surface tension. 

Experimentally, binary mixtures are prepared at various fixed compositions and the sur- 
face tension may then be observed as a function of temperature: see, e.g., [3]. This is 
comparable to keeping the order parameter M fixed. Accordingly we present such plots for 
the EdGF predictions in Fig. 

It must be recalled, however, that EdGF theory yields only the scaling or singular part of 
the surface tension. Since the leading power-law of the surface tension, vanishes at the 
critical point, the contribution from the analytic background is non-negligible. To illustrate 
this point, we include in Fig. Efb) plots for the full surface tension, E = AS + S , with an 
assumed but reasonably realistic model background term, namely, 

E {T,h)/K = 1 - 2t + 5t 2 , (6.13) 

where, for simplicity, only the temperature dependence has been considered. Since K + < 
[see (|6.6p ]. the scaling part of the surface tension above T c in zero-field, namely, &E a \py ~ 
K + \t\ >1 , curves downwards as seen in Fig. 0Ja). However, in Fig. Efb), the corresponding 
curvature in the full surface tension now appears to be upwards owing to the effects of the 
background So(T). In fact, the same sign of apparent curvature is observed in the NWW 
experiment on isobutyric acid and water |lj]; this clearly demonstrates the importance of 
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FIG. 5: Plots of (a) the singular part of the critical endpoint surface tension as a function of the 
reduced temperature t for fixed values of M and (b) with the addition of the model background 
term Eq(T, h) given in H6.13j) for the same values of M. 
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the background E (T, h). T 
will be discussed elsewhere 



le significance of also introducing integral powers of h in (j6.12j) 
22 1 - 



C. The Complete- Wetting Singularity 



It has been predicted by Cahn [26[ that a logarithmic singularity should occur in the 
slope of AS Q | 7 (T, h) on approaching the coexistence curve below T c , i.e., by taking h — > 0+ 
(along any generic, nontangential path). This logarithmic singularity is associated with the 
complete-wetting transition that, in turn, is reflected in Antonow's rule. Recall that when 
h = 0, the two liquid phases j3 and 7 coexist with the vapor phase a. Thus, when h — > 0+ 
while a and 7 phases coexist [i.e., on the surface a in Fig. [Ha)], the (3 phase of intermediate 
density emerges and spreads over (or wets) the a|7 interface. Indeed, one can exp licitly 
show that such a singularity arises within the EdGF theory: see Appendix C of |l7[. The 
analysis establishes that the coefficient of the ln/ioo term in (dAT,/dh) is positive, which 
is fully consistent with the numerical calculations presented in Fig. EJ these demonstrate 
a ln(T — To) -1 singularity in the derivative of the surface tension when T approaches the 
coexistence curve at fixed density p > p c , i.e., M > 0. 

Although the presence of a logarithmic singularity in £ a | 7 (/i) on approaching a wet in- 
terface may be regarded as well established theoretically, an important caveat is that in all 
cases the corresponding theoretical analysis entails the assumption that interactions within 
the fluids are entirely of a short-range character, i.e., decaying faster than any power- law. 
This, of course, precludes slowly decaying interaction potentials such as the 1/r 6 form that 
characterizes the van der Waals intermolecular forces prevalent in real molecular fluids. At 
fixed temperature T in the range TV < T < T c van der Waals forces should generate a 



(h - h c 



-1/3 



singularity in E a i 7 (/i). Thus, without special allowance for power-law poten- 



tials, a local functional analysis such as our EdGF theory must be suspect in relation to 
real fluid systems unless short range interactions happen to dominate for, say, reasons of 
symmetry, or accidental near-cancellation, etc. 

On the other hand, in the case of d = 3 bulk critical behavior it is known that 1/r 6 
potentials enter the asymptotic scaling forms only as irrelevant corrections-to-scaling (even 
though these may be dangerously irrelevant for certain quantities such as correlation func- 
tions at long distances) [471] . It is possible that a similar situation pertains in the case of the 
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FIG. 6: (a) Temperature derivative of AS = AS/K at fixed M (with values the same as in Fig. 03) 
revealing a logarithmic singularity when M is positive and t (< 0) approaches the coexistence curve; 
(b) semilogarithmic plot of the derivative of the scaling function S-(h) beneath T c : the straight 
line is a guide to the eyes while the dots represent the derivative. 
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surface tensions near a critical endpoint in which case a local functional theory might still 
prove asymptotically adequate. To our knowledge, however, this issue remains open and, as 
yet, a singularity in the derivative of the surface tension E Q | 7 (T, h) near a critical endpoint 
has not been identified experimentally: see, e.g., 



VII. SUMMARY 

As illustrated in Fig.Q binary fluid mixtures exhibit critical endpoints where, in the three- 
dimensional thermodynamic field space, a critical line of mixing transitions terminates at a 
first-order transition surface between the liquid phases and their common vapor phase. At 
such critical endpoints, the interfacial or surface tension S(T, h) becomes singular in both 
temperature and ordering field h. Our aim here has been to calculate the scaling functions 
describing the asymptotic behavior of the surface tensions through the whole neighborhood 
of the critical endpoint. 

To this end, the local functional theory of Fisher and Upton 0, 0] in the extended 
de Gennes-Fisher (EdGF) version of the theory, has been exploited because it captures 
many significant physical features tied to the nonclassical values of the relevant critical 
exponents, especially 77 > 0. As seen in f|4.1|) . the EdGF theory requires suitable scaling 
representations for the auxiliary free energy W(M; T, h^, g) and for the correlation length 
factor £ 2 /2x- To generate these, we have used the extended sine model of Ref. 9, since it 
embodies the appropriate analytic behavior, extends smoothly through the two-phase region 
below T c , and fits the values of many important universal amplitude ratios. 

However, both numerical and analytical calculations lead to the prediction of a small 
but unphysical cusp in the variation of the surface tension E(T, h) on crossing the critical 
isotherm at positive h, i.e., on entering the 7 region of the phase diagram: see Fig. ^a). 
This unanticipated behavior represents a shortcoming of the EdGF theory that is found to 
originate in the predicted variation of the order parameter profile, M(z), in the immediate 
vicinity of T = T c when it passes through the critical value M(zq) = 0. An improved 
local functional theory might avoid this difficulty To that end, it may be worthwhile to 
investigate the generalized or GdGF theory jjj]. However, the undesirable feature may still 
remain since we are inclined to believe that the origin of the cusp is associated rather directly 
with the single, scalar order-parameter formulation: both EdGF and GdGF theories employ 
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a simple scalar order parameter which cannot avoid "local criticality" when the interfacial 
profile crosses from one phase to another at T = T C . In response to this observation, Mikheev 
and Fisher 48] have addressed the formulation of two-order-parameter theories in which, 
in particular, the local energy fluctuation, as a second 'critical density,' plays a role; but a 
practicable scheme of approximation has not so far been achieved. 

Nevertheless, the application of EdGF theory to other properties of fluid interfaces and 
surfaces, seems worthwhile (e.g. juij]) and in the present case the cusp in s(8), the scaling 
function for the surface tension, produces a deviation of only a few percent from the naturally 
interpolated analytic variation. Accordingly, for numerical purposes we have adopted a 
smoothing procedure that removes the cusp; this yields the ameliorated angular function 
s(9) that is recorded numerically in Table I. 

On this basis the universal scaling functions S±(h) and S^(m) have been calculated: see 
Figs. |U and 03 which reveal significant features of the anticipated variation of the surface 
tension S(T, h). In particular, the role of the analytic background contribution So(T, h) 
can be assessed. Together with the improved theoretical predictions 1)6.8)1 for the universal 
surface tension amplitude ratios, P and Q, these results will be used elsewhere j^J to 
reanalyze the experimental data of Nagarajan, Webb, and Widom and to compare with 



other experiments 



18. 
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